Scattering theory of the hyperbolic BC n Sutherland and 
the rational BC n Ruijsenaars— Schneider— van Diejen 

models 



B.G. Pusztai 
Bolyai Institute, University of Szeged, 
Aradi vertaniik tere 1, H-6720 Szeged, Hungary 
'_^: ■ e-mail: gpusztai@math.u-szeged.hu 

o 

(N 

Abstract 

"^ ■ In this paper, we investigate the scattering properties of the hyperbolic BC n Suther- 

Q^ , land and the rational BC n Ruijsenaars-Schneider-van Diejen many-particle systems with 

three independent coupling constants. Utilizing the recently established action-angle du- 
ality between these classical integrable models, we construct their wave and scattering 
*~h, maps. In particular, we prove that for both particle systems the scattering map has a 
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1 Introduction 

The Calogero-Moser-Sutherland (CMS) many-particle systems (see e.g. [I], [2], [3], [I]) and 
their relativistic deformations, the Ruijsenaars-Schneider-van Diejen (RSvD) models (see e.g. 
[5], [6]) are among the most actively studied integrable systems. They appear in several branches 
of mathematics and physics, with numerous applications ranging from symplectic geometry, 
Lie theory and harmonic analysis to solid state physics and Yang-Mills theory. The intimate 
connection with the theory of soliton equations is a particularly important and appealing feature 
of these finite dimensional integrable systems. It is a remarkable fact that the CMS and the 
RSvD models associated with the A n root system can be used to describe the soliton interactions 
of certain integrable field theories defined on the whole real line (see e.g. [5], [7J, [8], |9], [TU]). 
In particular, these particle systems are characterized by conserved asymptotic momenta and 
factorized scattering maps (see e.g. [7], [11], [12], [13]). That is, in perfect analogy with the 
behavior of the solitons, the n-particle scattering is completely determined by the 2-particle 
processes. However, apart form some heuristic arguments [H] , the link between the integrable 
boundary field theories and the particle models associated with the non-A n -type root systems 
has not been developed yet. In our paper [T5] we have proved that the classical hyperbolic C n 
Sutherland model also has a factorized scattering map, but to our knowledge analogous results 
for the other non-A n -type root systems are not available. Motivated by this fact, in this paper 
we work out in detail the scattering theory of the hyperbolic BC n Sutherland and the rational 
BC n RSvD models with the maximal number of independent coupling constants. 
Upon introducing the subset 

c={i=(i 1 ,...,x n )6l n |i 1 >...>i n >0}cR n , (1.1) 

let us recall that the classical hyperbolic BC n Sutherland and the rational BC n RSvD models 
live on the phase spaces 

V s = {(q,p)\qe c, peW 1 } and V R = {(A, 6) \ A G c, 9 G M n }, (1.2) 

respectively. By endowing these spaces with the natural symplectic forms 

n n 

s - J2 d ^ A d p* and uR = Yl dXa A d ^' ( L3 ) 



a=l a=l 



we may think of V and V as two different copies of the cotangent bundle T*c. The hyperbolic 
BC n Sutherland model is characterized by the interacting many-body Hamiltonian 

n / 2 \ 

H S = J2[j + 9l^(Qa)+g 2 2 w(2q a )\+ ^ (g 2 w(q a - q b ) + g 2 w(q a + q b )) , (1.4) 

0=1 ^ ' l<a<b<n 

with potential function w(x) = sinh(x)~ 2 . To ensure the repulsive nature of the interaction, on 
the real coupling parameters g, g\ and g 2 we impose the constraints g 2 > and g\ + g\ > 0. 
As for the rational BC n RSvD model, the dynamics is governed by the Hamiltonian 



n n / A 2 \ 

i/« = ^oo S h(29„X(A) + ^n( 1 + ^) 

o=l ^ 0=1 V a J 



UK 

v 



;i.5) 



where the real parameters /x, v and k appear in the function 

too. Let us note that on the RSvD coupling parameters we impose the conditions \i ^ 0, v ^ 
and uk > 0. 

Working in a symplectic reduction framework, in our paper [16] we established the action- 
angle duality between the hyperbolic BC n Sutherland and the rational BC n RSvD models, 
provided the coupling parameters satisfy the relations 

9 2 = » 2 , #1 = 2^' 9l=^{v -nf. (1.7) 



By continuing our work [16], in this paper we explore the scattering theory of these particle 
systems. More precisely, after a brief overview in section [2] on the necessary background mate- 
rial, in section [3] we examine the temporal asymptotics of the trajectories of the Sutherland and 
the RSvD many-body models. The outcome of our analysis is formulated in lemmas [T] and [3J 
As expected, the asymptotics naturally lead to the wave and the scattering maps, too. Quite 
remarkably, the symplecticity of the wave maps can be seen as an immediate consequence of 
the duality between these integrable many-body models. We also prove that the scattering 
maps have factorized forms, which is the characteristic property of the soliton systems. Our 
main results on the wave and the scattering maps are summarized in theorems H] and HJ To 
conclude the paper, in section 0] we discuss our results and offer some open problems related to 
the scattering theory of the CMS and the RSvD systems. 

2 Preliminaries 

One of the main ingredients of the scattering theory we wish to develop in this paper is the 
action-angle duality between the hyperbolic BC n Sutherland and the rational BC n RSvD sys- 
tems [IB]. To keep the paper self-contained, in this section we gather some relevant facts about 
the symplectic reduction derivation of these models with emphasis on their duality properties. 
Up to some minor changes in the conventions, our overview closely follows [16]. For conve- 
nience, we formulate our results using the three independent RSvD coupling parameters /i, v 
and k. Without loss of generality, we may assume from the outset that v > > \x and k > 0. 

2.1 Background material from symplectic geometry 

Upon taking an arbitrary positive integer n e N, we introduce the set 

N„ = {l,...,n}cN. (2.1) 

In the following we also frequently use the notation iV = 2n. Now, with the aid of the unitary 
N x N matrix 

Un J-n 



(2-2) 



we define the non-compact real reductive matrix Lie group 

G = U{n, n) = {ye GL{N, C) | y*Cy = C} (2.3) 

with Lie algebra 

= u (u, n) = {Y e gl(N, C) | Y*C + CY = 0}. (2.4) 

Turning to V = G x g, at each point (y, Y) G V we define the 1-form 

# iy> Y)(Sy © SY) = tr (Yy^Sy) (Sy G T y G, 6Y G T y q = g). (2.5) 

It is clear that the product manifold V equipped with the symplectic form u = — d$ provides 
a convenient model for the cotangent bundle T*G. 

Next, for each column vector V G C^ we define the N x N matrix 

Z{V)=iii(W*-l N )+i(ii-v)C. (2.6) 

Also, by taking the unitary elements, in G (12.31) we choose the maximal compact subgroup 

K = U(n,n)f]U(N)^U(n)xU(n). (2.7) 

In its Lie algebra £ = u(n, n) fl u(N) the subset 

O = {£(V) \V eC N , V*V = N, CV + V = 0} c I (2.8) 

forms an orbit under the adjoint action of K. Consequently, it comes naturally equipped with 
the Kirillov-Kostant-Souriau symplectic structure uj° having the form 

u°([X,p],[Z,p])=tr(p[X,Z]) : (2.9) 

where p G O and [X, p] , [Z, p] G T p O with some X, Z G t. 

Our study of the CMS and the RSvD particle systems is based on the symplectic reduction 
of the extended phase space "P cxt = V x O endowed with the symplectic form 

u ext = ^(tu + tu°). (2.10) 

In passing we mention that the factor 1/2 is inserted into this definition purely for convenience. 
Now, let us note that the symplectic left action of the product Lie group K x K on V ext defined 
by the formula 

(k L ,k R ) ■ (y,Y,p) = (kLyk^^RYk^^Lpkl 1 ) ((k L ,k R ) G K x K) (2.11) 

admits a K x i^-equivariant momentum map. Making use of the identification (£ © t)* = £ © t 
induced by an appropriate multiple of the tr functional, the momentum map can be written as 

jext. pcxt _> | |? (y )Y) p) ^ ((yYy-% + p) © (-y f - wC), (2.12) 

where Y$ denotes the anti-Hermitian part of the matrix Y. 

The point of the above discussion is that by starting with certain K x .ft'-invariant Hamil- 
tonians on p^*, both the hyperbolic BC n Sutherland and the rational BC n RSvD models 
with three independent coupling constants can be derived by reducing the symplectic manifold 
(P ext ,u; cxt ) at the zero value of the momentum map J ext . 



2.2 The Sutherland model from Marsden— Weinstein reduction 

One of the key objects in the symplectic reduction derivation of the hyperbolic BC n Sutherland 
model is its Lax matrix L : V s — > g. It has the form 



A B 
-B -A 



-IkC, (2.13) 



where A and B are appropriate n x n matrices. More precisely, their entries are 

K, b = -* , , A v = Pa B a , b = —St- Bv = | " + "°°eW , (2.14) 

smh(g a - q b ) smh(g a + q b ) smh(2g c ) 

where a, b, c G N n and a ^ b. Next, to each real n-tuple q — (qi, . . . , q n ) we associate the matrix 

Q = diag(gi, . . . , q n , -q u . . . , -g n ), (2.15) 

and we also introduce the column vector E e C^ with components 

E a = -E n+a = 1 {ae N„). (2.16) 

Finally, let ?7(1)* denote the diagonal embedding of U(l) into K x K and define the product 
manifold 

M s = P s x(ifxif)/[/(l)». (2.17) 

Having equipped with the above objects, we can start the reduction procedure by analyzing 
the level set 

£ = {{y,Y,p) e V cxt | J cxt (y, Y,p) = 0}. (2.18) 

It turns out to be an embedded submanifold of V ext , and the diffeomorphism 

T S :M S ^ £ , (q,P, (vl,Vr)U(1)*) H- (^e^^ 1 ,^L(g,p)^ 1 ,^e(^)% 1 ) (2.19) 

gives rise to the identification .M 5 = -Co- By inspecting the (residual) K x K- action (12. lip on 
M. s , it is clear that the base manifold of the trivial principal (K x K)/U(l) ^-bundle 

7t s :M S ^ V s , (q,p, (til,VrM1)*) H- (?,p) (2.20) 

provides a realization of the reduced manifold. Since (tt s )*u s = (Tp)*ci; ext , we conclude that the 
Sutherland phase space V s (11.2)1 does serve as a convenient model for the reduced symplectic 
manifold V ext //o(K x K). Moreover, making use of the K x K- invariant Hamiltonian 

F: V cxt -+ R, (y,F,p) ^ F(j/,y,p) = ^tr(r 2 ), (2.21) 

we find the relationship (n s )*H s = (Tq)*E as well. In other words, the reduced Hamiltonian 
induced by the quadratic function F coincides with the Hamiltonian (ll.4p of the hyperbolic 
BC n Sutherland model with coupling constants displayed in equation (II. 7p . 



2.3 The symplectic reduction derivation of the RSvD model 

Compared to the Sutherland model, the reduction derivation of the rational BC n RSvD model 
requires a slightly longer preparation. As a first step, for each a G N n we define the function 

, a ^^,-( 1+ gn(i +3 ^)(i +3 ^),c „ 22) 

Also, we introduce the function A: V R — > G with matrix entries 



A p "a "b I y y I o A p 



2i/i „ _ 0a+0t z a z 6 2i/i 



2i/x + A a - A 6 ' * "-■-'•'■"-'-' |z a z 6 |2 2i/i - A a + A 



b 



A -~A -r-Oa+Ot-U --1|* 2i ^ | K^-^ x 



(2.23) 



where a,b G N n . As we have shown in our earlier paper [T7], the positive definite N x N matrix 
A provides a Lax matrix for the rational C n RSvD model with parameters // and v. 
Next, with the aid of the K-dependent functions 



a(x) = V x + Vx 2 + k 2 and g(a . ) = i/6 1 1 (2>24) 

V2x V2^ y/x + v^ 2 + /t 2 



G G. (2.25) 



defined for x > 0, for each A G c we introduce the N x N matrix 



h{\) 



diag(a(Ai), . . . , a(X n )) diag(/3(Ai), . . . , /3(A n )) 
-diag(/3(Ai), . . . , /3(A n )) diag(a(Ai), . . . , a(A n )) 



Utilizing /i(A), the Lax matrix of the rational BC n RSvD model can be written as 

A BC : V R -)• G, (A, 0) H- „4 BC (A, 0) = /i(A)-U(A, 9)h(\)-\ (2.26) 

Heading toward our goal, we still need the column vector J-"(A, #) G C^ with components 



T a (\,9) = e-°«\z a (\)\* and F n+a (\,6) = e y °z a (\)\z a (\)\-^ (2.27) 

where a G N n , together with the column vector 

V(\,6)=A(\,6)-^(\,9)eC N . (2.28) 

Also, for each A = (Ai, . . . , A n ) G M. n we define the diagonal matrix 

A = diag(Ai, . . . , A n , -A 1; . . . , -A n ). (2.29) 

At this point we can put the RSvD particle system into the context of symplectic reduction. 
Upon introducing the smooth product manifold 

M R = V R x (K x K)/U(l)*, (2.30) 



the diffeomorphism T^: Ai R — >■ £o defined by the formula 

(A, 9, ( Vl , 7i R )U(l)*) ^ ( V lA(X, 9)h(X)- 1 rj^\ VR h(X)Ah(X)- 1 ^ 1 , 7^(V(A, O))^ 1 ) (2.31) 

provides an alternative parametrization of the level set £o A2.18J) . Moreover, the explicit form 
of the (residual) K x fC-action ( 12 .lip on the model space M. R = £ permits us to identify the 
reduced manifold with the base manifold of the trivial principal (K x K)/U(l)*-bwa.dle 

n R : M R _» p* ( A , 9, (r] L , r} R )U{l) # ) H- (A, 9). (2.32) 

Since (n R )*uj R = (T^)*o; ext , we conclude that the reduced space V ext //o(KxK) can be naturally 
identified with the RSvD phase space V R (II. 2p . Finally, let us consider the function 

/ : p-t ^ Rj (y5F5p) ^ /(^y^) = ltr(yy*). (2.33) 

From the relationship (7r R )*H R = (T R )*f we infer immediately that the reduced Hamiltonian 
function corresponding to the K x f^-invariant function / coincides with the Hamiltonian ( II. 5p 
of the rational BC n RSvD model. 

2.4 Action-angle duality 

The two different parameterizations Tq ( I2.19P and T^ (I2.3ip of the level set £o (I2.18P induce 
two different models, V s and V R fll.2[) . of the symplectic quotient V ext //o{K x K). Therefore 
there is a symplectomorphism 

S:V S ^V R (2.34) 

between V s and V R , uniquely characterized by the equation 

S o n s o (T^)- 1 = 7i R o (T^)" 1 . (2.35) 

Making use of this purely geometric observation we can easily construct canonical action-angle 
variables for both the Sutherland and the RSvD models. Indeed, by pulling back the canonical 
positions and momenta of the RSvD system, the global coordinates S*X a and S*9 a on V s give 
rise to canonical action-angle variables for the Sutherland model. Similarly, by pulling back the 
canonical positions and momenta of the Sutherland system, the global coordinates (iS _1 )*g a and 
(S~ 1 )*p a on V R provide canonical action-angle variables for the RSvD model. This remarkable 
phenomenon goes under the name of action-angle duality between the Sutherland and the RSvD 
particle systems. 

3 Scattering theory 

In this section we provide a thorough analysis of the time evolution of the hyperbolic BC n 
Sutherland and the rational BC n RSvD dynamics for the large positive and negative values of 
time t. Based on the resulting temporal asymptotics we construct the wave and the scattering 
maps as well. 
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3.1 Scattering properties of the Sutherland model 

Take an arbitrary point (q, p) G V s and consider the Hamiltonian flow of the Sutherland model 

R3t^(q(t),p(t))EV S (3.1) 

satisfying (q(0),p(0)) = (q,p). Recalling the Sutherland Lax matrix ( I2.13p . we let L = L(q,p). 
Keeping the notation (" 12. 15[) in effect, for each teRwe also define 

Q(t) = diag( 9l (t), . . . , q n (t), - qi (t), . . . , -q n (t)). (3.2) 

Now, it is clear that the (complete) 'geodesic' flow 

R9t^(e Q e fL ,I,((£))e£ (3.3) 

generated by the unreduced free Hamiltonian F (12.2 ip projects onto the (complete) reduced 
flow ( 13. ip . In particular, utilizing the parametrization Tq ( I2.19p . for all t G M we can write 



QJ.L 



e^e 



k^e^k^t)- 1 (3.4) 



with some ki(t), knit) G K. It entails the spectral identification 

a(e 2Q ^) = a(e 2Q e tL e tL *) } (3.5) 

which can be seen as the starting point of a purely algebraic solution algorithm of the Sutherland 
model based on matrix diagonalization. 

Making use of the symplectomorphism S (I2.34p , the above matrix flows can be parametrized 
with the dual variables (A, 6) = S(q,p) as well. First, recalling the matrices (I2.23J) and ( I2.25p . we 
introduce the shorthand notations A = ^.(A,^) and h = h(X). Second, using the abbreviation 
defined in (12.291) . notice that the parametrization T^ (I2.3ip together with the defining property 
of S (12.351) immediately lead to the relationships 

e Q =i] L A^h~ 1 ri R 1 and L = rjRhAh' 1 ^ 1 (3.6) 

with some tjljTJr £ K- ft readily follows that 

e 2Q e tL e tL * = m h~ x Ae tk h~ 2 e tK hif R \ (3.7) 

However, since the functions (I2.24p appearing in the definition of h ( I2.25P obey the identities 



K 2 IK 



a(xY + /3(xY = 1, a(xY - f3(xY = \ 1 + — , 2a(x)/3(x) = — , (3.8) 

V x 1 x 

we find easily that 

h- 2 = y/l N + k 2 A- 2 + ikCA- 1 . (3.9) 

Therefore, by plugging the above formulae into (13.51) . we obtain 

a(e 2Qit) ) = a [a^1 n + K 2 A~ 2 e 2tA + inACA- 1 ) . (3.10) 



By exploiting the consequences of this formula, now we can work out the scattering theory of 
the Sutherland model. Before going into the details, for each a G N n and A G c we define 



^>HiH+(^)<4K-(^y 



+5EM 1 



6=1 



x b y 
V 



1 / V Z \ 1 / K Z 



(3.11) 



Lemma 1. Take an arbitrary (q,p) G V s and let (X,9) = S(q,p) G V R . Consider the Hamil- 
tonian flow of the Sutherland model 



3t^(q(t),p(t))eV' 



(3.12) 



satisfying (q(0) , p(0)) = (q,p). Then there are some positive constants T > 0, r > and C > 0, 
such that for all a G N n and for all t > T we have 



q a (t) -(t\ a -9 a + -A a (X) 



< Ce- tr and \p a (t) - X a \ < Ce 



-tr 



whereas for all a G N ra and for all t < —T we can write 



q a (t)-(t(-X a )+9 a + -A a (X) 



< Ce tr and \p a (t) - (-X a )\ < Ce 



tr 



(3.13) 



(3.14) 



Proof. To prove the lemma, our guiding principle is Ruijsenaars' result on the temporal asymp- 
totics of the spectra of exponential matrix flows (see Theorem A2 in [7]). First, we introduce 
the n x n matrix TZ n with entries 

(J^n)a,b = 8a+b,n+l- (3.15) 

By conjugating with 



w 



1Z 



GGL(AT,C), 



(3.16) 



we also define the N x N matrices 

M = WA^1 N + k 2 A- 2 W _1 , D = WAW" 1 , X = ikW AC ^ l W~ l . (3.17) 

Let us observe that the diagonal matrix D has the property that its eigenvalues are in strictly 
decreasing order on the diagonal. 

By inspecting the matrix entries of M (I3.17p . notice that for all a, b G N n we have 



M a , 6 = Aj(l + k 1 X 



2\-2\ 



(3.18) 



For each a G N n let M^ denote the leading principal a x a submatrix taken from the upper- 
left-hand corner of M. Since M^ is essentially a Cauchy matrix, for its determinant we have 



det(M^) = l[e- 2e »\z b (X)\(l + K 2 X^ }[ (1 + 4/i 2 (A c - A c 



-2\-l 



(3.19) 



6=1 



l<c<d<a 



By taking the quotients of the consecutive leading principal minors of M, we also define 

det(M^) 
m 1 = M x x and m a = - ,\. ' (2 < a < n). (3.20) 

det(M(«- 1 )) \ - - ) \ ) 

From equation (13.191) it follows that for all a £ N n we can write 

a-l 

e- 2 ^|z a (A)|(l + k 2 \- 2 )1 \[{l + V(A a - A,)- 2 )- 1 . (3.21) 



a-l 
— 9fl- i l a \ I /* . 9 i _9s I 

6=1 



Therefore, recalling the formulae (12.22)) and (13.1 ip . we end up with the concise expression 

ln(m a ) = -2# a + A a (A). (3.22) 

Utilizing the above objects, we are now in a position to analyze the asymptotic properties 
of the trajectory (13.121) . By combining equations (13.101) and (13. 17)) . it is clear that for all t £ M 
we can write 

a(e 2Q ^) = a(Me 2tD + X). (3.23) 

By slightly generalizing Ruijsenaars' aforementioned theorem, one can easily verify that the 
exponentially growing large eigenvalues of the matrices Me 2tD and Me 2tD + X have essentially 
the same asymptotic properties as t — > oo. More precisely, there are some positive constants 
T > 0, r > and C > 0, such that for all a £ N n and for all t > T we have 

e 2 ^=m a e 2tx «(l + e a (t)) } (3.24) 

where the error term e a (t) obeys the estimation 

max{|e a (t)|,|e a (t)|} <Ce~ tr <\. (3.25) 

By taking the logarithm of equation (13.24)1 . it readily follows that 



Q a (t) -t\ a --\n(m a ) 



^\Hl + e a {t))\<\e a {t)\<Ce- tr . (3.26) 



Due to the formula (13.221) . for the large positive values of t the control over q a (t) is complete. 
Next, by taking the derivative of equation (13.241) . we also find 

However, the Hamiltonian equations of motion yield p = q, whence we obtain 

\p(t) - AJ = - , l £a ^J < |e a (*)| < Ce~ tr . (3.28) 

IfaV J a| 2 |l + e a (t)| ~ ~ 
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To conclude the proof, let us observe that by conjugating with the N x N matrix H N (j3.15p . 
the spectral identification (I3.23P can be rewritten as 

a(e 2Qit) ) = a {lZ N M1Z N l e 2tnNDn ^ + K N XKjf\ . (3.29) 

Since the eigenvalues of the diagonal matrix IZ^DIZ^ = —D are in strictly increasing order on 
the diagonal, the asymptotic relationships f l3.14|) for t — > — oo can be established by the same 
techniques that we used for the t — > oo case. □ 

Keeping the notations introduced in lemma [TJ the asymptotics can be rewritten as 

Q a (t)~tpt + qt and Pa (t)~Pa (<->±oo), (3.30) 

where the asymptotic phases and momenta have the form 

<£ = =F0a + ^A a (A) and p± = ±\ a , (3.31) 

respectively. Notice that the asymptotic states (q ± ,p ± ) belong to the manifolds 

V ± = {{x, y) = (x u ...,x n ,y 1 ,...,y n )eR 2n \y l ^...^y n ^ 0}, (3.32) 

that we endow with the natural symplectic forms 



w* 



^2 dx a A dy a . (3.33) 



o=l 



One of the principal goals of scattering theory is to study the wave maps 

Wl:V S -> V ± , (q,p) H- (q ± ,p ± ). (3.34) 

Theorem 2. The wave maps W± of the Sutherland model are symplectomorphisms from V s 
onto V^ . The scattering map S s = Wf o (W_) _1 is also a symplectomorphism of the form 

V- 3 (x, y) i-> S s (x, y) = {-x x + Ai(-y), ...,-x n + A n (-y), -y x , . . . , -y n ) E V + . (3.35) 

Proof. Let us introduce the maps 

T s . v r ^ p±) (A; e) ^ f T0i + I Al(A ) 7 . . . ? T0n + lA n (A), ±Ai, . . . , ±X n ) . (3.36) 



Recalling A a (13. lip , it is evident that 7~± are symplectomorphisms with inverses 

(TlY\x, y) = ( ± Vl , ..., ±y n , txi ± \&i{±y), ■■-,Tx n ± l&n(±y) J • (3.37) 



Moreover, remembering the asymptotic phases and momenta (13.31 j) . it readily follows that the 
wave maps ( I3.34p are symplectomorphisms of the form 

W£ = T£ o S. (3.38) 

Since S s = T+ o (7^) _1 , the explicit formula ( I3.35P is also immediate. D 
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3.2 Scattering properties of the RSvD model 

Take an arbitrary (A, 9) G V R and consider the Hamiltonian flow of the RSvD model 

R3t^ (X(t),0(t)) eV R (3.39) 

passing through the point (A, 9) at t — 0. Remembering the matrices (I2.23p . (12.25 p . (I2.26p . and 
the column vector (I2.28p . we introduce the abbreviations 

A = A(\,9), h = h(\), A BC = A BC (\,9), V = V(X,9). (3.40) 

Keeping the notation displayed in equation (12. 29 p . for all £ G M we also the define diagonal 
matrix 

A(£) = diag (Ai(t), . . . , A„(£), -Ai(t), . . • , -K{t)) ■ (3.41) 

Now, bearing in mind the symplectic reduction derivation of the RSvD model, let us note that 
the (complete) 'linear' flow 

R 3 t H- (Ah' 1 , hAh- 1 - t (A BC - (A 30 )' 1 ) , £(V)) G £ (3.42) 

generated by the unreduced Hamiltonian / (I2.33P projects onto the (complete) reduced RSvD 
flow (I3.39p . Recalling the parametrization T R (I2.3ip . it is evident that for all t & R we have 

hAh' 1 - t (A BC - (A 30 )- 1 ) = A; iJ (t)/ i (A(t))A(t)/i(A(t))- 1 A; K (£)- 1 (3.43) 

with some kn(t) G K. In particular, we obtain the spectral identification 

a(A(t)) = a (hAh' 1 - 1 (A BC - (A 30 )- 1 )) . (3.44) 

The point is that each trajectory of the RSvD dynamics can be recovered by simply diagonal- 
izing a linear matrix flow. 

Utilizing the duality symplectomorphism S f )2.34p . our next goal is to parametrize the above 
linear matrix flow with the dual variables (q,p) = 5 _1 (A, 9) G V s . Recalling the Sutherland Lax 
matrix (I2.13p . we first define L = L(q,p). At this point, remembering the notation (I2.15p . it is 
clear that the form of Tq ( I2.19P and the defining property of S ( I2.35P lead to the relationships 

A^h' 1 = ri L e Q rf^ and hAh" 1 = tjrLt]^ 1 (3.45) 

with some r]L,r}R G K. It entails A BC = ^Re 2 ^^ 1 , therefore the matrix 

^bc _ (^BC)-i = 2r]R sinh ( 2 Q)^i (3.46) 

has a simple spectrum. Moreover, recalling the relationship (I3.44p . for the spectrum of A(t) we 
obtain the particularly useful formula 

a(A(t)) = a(L-2t sinh(2Q)) . (3.47) 
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Lemma 3. Take an arbitrary point (X,9) G V R and let (q,p) = S 1 (A,^) G V s . Consider the 
Hamiltonian flow of the RSvD model 

R3t^(\(t),6(t))eV R (3.48) 

satisfying (A(0), #(0)) = (A, 6). Then there are some positive constants T > and C > 0, such 
that for all a G N ra and for all t > T we have 

\X a (t) - (2tsmh(2q a )-p a )\ < Cr 1 and \O a {t) - q a \ < Ct~ 2 , (3.49) 

whereas for all a G N„ and for all t < —T we can write 

\X a (t)-(2tsmh(-2q a )+p a )\<C\t\- 1 and \0 a (t) - (-q a )\ < C\t\- 2 . (3.50) 

Proof. Making use of lZ n ( I3.15P , we define 

GGL(A,C). (3.51) 



W 



On lr 

n n o r 



With the aid of the N x N matrices 

D = -2W sinh(2Q)W- 1 and M = WLW~\ (3.52) 

the spectral identification (I3.47P can be cast into the form 

<r(A(t)) = (r(£D + M), (3.53) 

where D is a diagonal matrix with its eigenvalues in strictly decreasing order on the diagonal. 
Therefore elementary perturbation theory can be readily applied to analyze the properties of 
A(t) as t — > oo. (For a short account on the relevant facts from perturbation theory see e.g. 
Theorem Al in [7].) More precisely, there are some positive constants To > and C > 0, such 
that for all a G N n and for all t > T the eigenvalue A a (t) has the form 

X a (t) = tD aA + M a>a + e a (t) = 2t sinh(2g a ) - p a + e a (t), (3.54) 

where the error term e a (t) satisfies 

|ea(*)|<Ct -1 and \i a {t)\ < Ct~ 2 . (3.55) 

Since the asymptotic property of X(t) is under control, now we can turn our attention to 
the asymptotic analysis of 0(t). Making use of the equations of motion generated by the RSvD 
Hamiltonian H R (jl.5p . we find 

X a {t) = 2 smh(20 a (t))v a (X(t)). (3.56) 

From the asymptotic behavior of X(t) f)3.54p it is clear that there are some constants 7\ > To 
and K, > 0, such that for all a, 6, c G N n , a ^ b, and for all t > T\ we have 

\X a (t) - \ b (t)\ > Kt and X c {t) > Kt. (3.57) 
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By inspecting v a ( \1.6\i . it is also evident that there are some constants T > T\ and H > 0, such 
that for all a£N„ and for all t > T we can write 

v a (X(t)) < l + Ur 2 . (3.58) 

By combining equations H3.54J) and (I3.56p . for all o£N„ and for all t > T we obtain 

• uoa (A\ ■ Uo \ (l-^(A(t)))sinh(2g a ) + |e a (t) 
smh(20 a (t)) - smh(2g a ) = — — — * , (3.59) 

Va{A-{t)) 



whence the estimation 



sinh(20 a (f))-sinh(2g a )| /^sinh(2gi) C\ 1 
"2 - V 2 + 4 j ^ 



|*.(t) - q a \ < l^^-^^l < ( "»"^i; + - ) ± 2 ( 3. 6 o) 



is immediate. Thus the asymptotic relationships (I3.49P are established. 

To conclude, let us note that by applying the same techniques on the matrix flow 

t ^ TZ N {tD + M)TZ]^ , (3.61) 

one can easily prove the remaining relationships (I3.5(jp as well. □ 

Keeping the notations of lemma [3J our results on the asymptotics can be rephrased as 

\ a (t) ~ 2tsinh(26t) + \t and 6 a {t) ~ 0j (t-»±oo), (3.62) 

where for all a G N n we have 

A± = T p a and 6£ = ±q a . (3.63) 

Notice that the asymptotic states (X ± , 6 I± ) belong to the phase spaces V ± (I3.32p . therefore the 
RSvD model can be characterized by the wave maps 

Wi:V R ^V ± , (X,9)^(X ± ,9 ± ). (3.64) 

Theorem 4. The wave maps W R of the RSvD model are symplectomorphisms from V R onto 
V ± . The scattering map is also a symplectomorphism of the form 

S R = W* o (W R Y l : V~ -)• V + , (z, y) ^ S R {x, y) = {-x, -y). (3.65) 

Proof. With the aid of the symplectomorphisms 

T£ : V s -» P ± , (?,p) H> (=FP, ±<?) (3.66) 

the wave maps (I3.64p can be realized as compositions of symplectomorphisms of the form 

W R = r R oS-\ (3.67) 

Due to the relationship S R = T R o (TJ^) -1 , the formula (I3.65P also follows. D 
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4 Discussion 

For many physical systems it is the scattering theory that provides the main tool to investigate 
the properties of the constituent particles and the nature of their interaction. At the same 
time, scattering theory is a notoriously difficult subject heavily relying on non-trivial techniques 
from hard mathematical analysis. However, by exploiting the action-angle duality between the 
Sutherland and the RSvD models, a careful analysis of their algebraic solution algorithms 
led us to rigorous temporal asymptotics of the trajectories. Having a glance at the structure 
of the resulting wave maps W± (I3.38J) and W± (13.671) . it is transparent that they are made 
of two strikingly different building blocks. The explicitly defined maps T± (13.361) and T± 
(13.661) are manifestly symplectic, meanwhile the construction and the symplecticity of S (I2.34p 
hinges on the symplectic reduction derivation of the Sutherland and the RSvD systems. That 
is, beside providing canonical action-angle variables, the geometric ideas culminating in the 
duality symplectomorphism S H2.34J) prove to be fundamental in the scattering theory as well. 
Looking at the scattering maps S s (13.351) and S R (I3.65p . we see that, up to an overall sign, 
the asymptotic momenta of both the hyperbolic BC n Sutherland and the rational BC n RSvD 
models are preserved. Following Ruijsenaars' terminology [13], we can say that these integrable 
systems are pure soliton systems of type BC n . Moreover, for the Sutherland model the classical 
phase shifts are completely determined by the 2-particle processes and the 1-particle scattering 
on the external field. As for the rational RSvD system, the phase shifts are trivial. In other 
words, in both cases the scattering map has a factorized form. Due to the puzzling connection 
with the theory of solitons, it would be of considerable interest to develop an analogous theory 
at the quantum level, too. Relatedly, we find it a much more ambitious, but equally motivated 
research problem to classify the pure soliton systems associated with arbitrary root systems. 
We wish to come back to these issues in later publications. 
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